Chaos and Chaotic Phase Mixing in Galaxy 
Evolution and Charged Particle Beams 



Henry E. Kandrup 

University of Florida, Gainesville, FL 32611, USA 

Abstract. This paper discusses three new issues that necessarily arise in realistic 
attempts to apply nonlinear dynamics to galaxy evolution, namely: (i) the meaning 
of chaos in many-body systems, (ii) the time-dependence of the bulk potential, which 
O ■ can trigger intervals of transient chaos, and (iii) the self-consistent nature of any bulk 

chaos, which is generated by the bodies themselves, rather than imposed externally. 
Simulations and theory both suggest strongly that the physical processes associated 
with galactic evolution should also act in nonneutral plasmas and charged particle 
beams. This in turn suggests the possibility of testing this physics in real laboratory 
experiments, an undertaking currently underway. 
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^ 1 Introduction 

o ' 

■ As recently as 1990, most galactic dynamicists ignored completely the possible 

role of chaos in galaxies. However, the past decade has seen a growing recognition 
C^l ■ that chaos can be important in determining the structure of real galaxies. Still, 

much recent work involving chaos in galactic astronomy has been simplistic in 
that it has involved naive applications of standard techniques from nonlinear dy- 
namics developed to analyse two- and three-degree-of-freedom time-independent 
q ■ Hamiltonian systems. The object here is to discuss some of the additional com- 

S—{ | plications which arise if nonlinear dynamics is to be applied to real galaxies, 

which are many-body systems comprised of a large number of interacting masses 
and characterised by a self- consistently determined bulk potential which, during 
their most interesting phases, can be strongly time- dependent. 
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H ' 2 Transient chaos and collisionless relaxation 

2.1 Transient chaos induced by parametric resonance 

It is well known to nonlinear dynamicists that the introduction of an oscilla- 
tory time-dependence into even an otherwise integrable potential can trigger an 
interval of transient chaos, during which many orbits exhibit an exponentially 
sensitive dependence on initial conditions. Physically, this transient chaos arises 
from a resonance overlap between the frequencies ~ Q of the unperturbed orbits 
and the frequency or frequencies ~ w of the time-dependent perturbation. 

In the past, the possible effects of such chaos have been considered for both 
nonneutral plasmas [jlj and charged particle beams |2j . More recently, such tran- 
sient chaos has also begun to be considered in the context of galactic astron- 
omy H. That work has shown that, for large fractional amplitudes, > 0.1 or 
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so, this resonance can be very broad, triggering significant amounts of chaos for 
10 _1 <u)/Q< 10; and that the existence of the phenomenon is robust, compar- 
atively insensitive to details. It will, for example, persist if one allows for damped 
oscillations and/or modest drifts in frequency (generated, e.g., by making u> a 
random variable sampling an Ornstein-Uhlenbeck colored noise process). 

The breadth of the resonance and the insensitivity to details suggest that 
transient chaos could well prove common, if not generic, in violent relaxation Q, 
the collective process whereby a (nearly) collisionless galaxy or galactic halo 
evolves towards an equilibrium or near-equilibrium state. Violent relaxation typ- 
ically involves damped oscillations triggered, e.g., by interactions with another 
galaxy or, in the early Universe, by the cosmological details preceding galaxy 
formation. However, when considering collective effects there is only one natural 
time scale, the dynamical time tp ~ 1/ ' \JGp, with p a typical mass density, which 
determines both the characteristic orbital time scale and (at least initially) the 
oscillation time scale. The exact numerical values of these time scales will involve 
numerical coefficients which will in general be unequal and vary as a function 
of location within the galaxy. If, however, one need only demand that the oscil- 
lation and orbital time scales agree to within an order of magnitude, it would 
seem likely that this resonance could trigger transient chaos through large parts 
of the galaxy. In real galaxies the oscillations will presumably damp and the 
frequencies drift as the density changes and, presumably, power cascades from 
longer to shorter scales. To the extent, however, that the details are unimportant 
such variations should not obviate the basic effect. 



2.2 Chaotic phase mixing and collisionless relaxation 

But why might such transient chaos prove important in galactic evolution? De- 
tailed numerical simulations indicate that violent relaxation can be a very rapid 
and efficient process, but simple models involving regular orbits, such as Lynden- 
Bell's H balls rolling in a pig-trough, do not approach an equilibrium nearly fast 
enough ||. The important point, however, is that allowing for the effects of chaos 
can in principle dramatically accelerate both the speed and efficacy of violent 
relaxation. An initially localised ensemble of regular orbits evolved into the fu- 
ture in a time-independent potential will begin by diverging as a power law in 
time and, only after a very long period, slowly evolve towards a time-averaged 
equilibrium, i.e., a uniform population of the K AM ton to which it is restricted. 
By contrast, a corresponding ensemble of chaotic orbits will begin by diverging 
exponentially at a rate that is comparable to a typical value of the largest fi- 
nite time Lyapunov exponent for the orbits in the ensemble; and then converge 
exponentially towards an equilibrium or near-equilibrium state at a somewhat 
smaller, but still comparable, rate. The exponential character of this chaotic 
phase mixing means that the time scale associated with this process is typically 
far shorter than the time scale associated with regular phase mixing [^|J7|J^] . 

It is evident that chaotic phase mixing in a time-independent Hamiltonian 
system can trigger a very rapid approach towards an equilibrium, but this does 
not necessarily 'explain' violent relaxation. If, e.g., most of the orbits in the 
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system are regular, it would seem unlikely that chaotic phase mixing could be 
sufficiently ubiquitous to explain an approach towards a (near-)equilibrium for 
the galaxy as a whole, fndeed, one would expect that, for a galaxy that is in or 
near equilibrium the relative measure of chaotic orbits should be comparatively 
small: If the galaxy exhibits nontrivial structures like a bar or a cusp, the types 
of structures which one has come to associate with chaos, one would also expect 
large numbers of regular orbits must be present to serve as a 'skeleton' to support 
that structure ||. Moreover, it is evident that, although chaotic mixing in a 
time-independent potential can be very efficient in mixing orbits on a constant 
energy surface, the energy of each particle remains conserved, so that there can 
be no mixing in energies. The extent to which chaotic phase mixing in a time- 
dependent potential will trigger an efficient shuffling of energies is not completely 
clear. 

The important point, then, is that chaotic phase mixing associated with 
transient chaos in a time-dependent potential is likely to explain these remain- 
ing lacunae. At least for large amplitude perturbations, (say) 10% or more, this 
parametric resonance can trigger a huge increase in the relative abundance of 
chaotic orbits so that, for pulsation frequencies near the middle of the reso- 
nance, virtually all the orbits exhibit substantial exponential sensitivity. More- 
over, given that this chaos involves a resonant coupling, it tends typically to 
cause a substantial shuffling of energies: those frequencies which arc most apt to 
trigger lots of chaos are also apt to induce the largest shuffling of energies. 
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Fig. 1. (top) The emittance e x for an ensemble of initial conditions evolved in an 
integrable Plummer potential subjected to damped oscillations with four different fre- 
quencies, (bottom) x-y scatter plots corresponding to the uppermost curve. 
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Still it should be noted that one can get a 'near-complete' shuffling of orbits 
on different constant energy surfaces even if the orbital energies are not especially 
well shuffled. This, however, is not necessarily a problem. Simulations of systems 
exhibiting efficient collective relaxation do not necessarily involve masses which 
completely 'forget' their initial conditions. Rather, comparatively efficient and 
complete violent relaxation is completely consistent with an evolution in which 
masses 'remember' (at least partially) their initial binding energies, i.e., in which 
masses that start with comparatively large (small) binding energies end up with 
comparatively large (small) binding energies JhJ. 

That it may be possible to achieve efficient chaotic phase mixing in an os- 
cillating galactic potential while still relaxing towards a nearly integrable state 
within lOtp or so is illustrated in Fig. |l|. This Figure was generated from orbits 
evolved in a time-dependent potential of the form 

m(t) . snout 

F( ^' M) = - (i + ^ + ^ + z^/2 ' mW = 1 + 5m (^T^' (1) 

with 5m — 0.5, t a — 100 and p = 2, which represents a galaxy damping towards 
an integrable Plummer sphere. The four curves in the top panel exhibit the x- 
component of the phase space emittance, = ((rf)(vf) — (r^Wi) 2 ) 1 / 2 (i = x, y, z), 
all computed for the same localised ensemble of initial conditions, but allowing 
for four different frequencies u>. The curves exhibit considerable structure but, at 
least for early times, the overall evolution is exponential. The bottom panels ex- 
hibit the x and y coordinates at five different times for the ensemble represented 
by the uppermost of the four curves. Here tjy ~ 20, so that t — 256 corresponds 
to roughly 12tn- 

Intuitively, one might expect that strong oscillations, which trigger the largest 
finite time Lyapunov exponents and the largest number of chaotic orbits, would 
yield the fastest chaotic phase mixing and, hence, the most rapid and most 
complete violent relaxation. A time-dependence with a weaker oscillatory com- 
ponent, e.g., a time-dependence corresponding to a near-homologous collapse, 
might instead be expected to yield less chaos and, hence, less efficient and less 
complete violent relaxation. There is, therefore, an important need to determine 
the extent to which, in real simulations of violent relaxation, many/most of the 
orbits (or phase elements) are strongly chaotic, and the degree to which the rate 
and completeness of the observed violent relaxation correlate with the size of the 
largest finite time Lyapunov exponents and/or the relative measure of chaotic 
orbits. Investigations of these issues are currently underway. 



3 The role of discreteness effects 
3.1 Microchaos and macrochaos 

The discussion in the preceding section, like most applications of nonlinear dy- 
namics to galactic astronomy, neglects completely discreteness effects associated 
with the 'true' many-body potential, assuming that masses in a galaxy can be 



Chaos in Galaxy Evolution and Charged Particle Beams 



approximated as evolving in a smooth, albeit time-dependent, three-dimensional 
potential and that 'chaos' has its usual meaning. That this is justified is not com- 
pletely obvious. The gravitational TV-body problem for a large number of bodies 
of comparable mass is strongly chaotic in the sense that individual orbits have 
large positive Lyapunov exponent \n even when there is absolutely no chaos 
in the continuum limit! If, e.g., a smooth density distribution corresponding to 
an integrable potential is sampled to generate an TV-body density distribution, 
one finds that orbits evolved in this TV-body distribution will be strongly chaotic, 
even for very large TV, despite the fact that characteristics in the smooth potential 
generated from the same initial condition are completely integrable. Moreover, 
there is no sense in which the exponential sensitivity decreases with increasing 
TV: if anything \n is an increasing function of TV jy]. In this sense, larger TV 
implies more chaos, not less! 

This situation has led some astrophysicists to question, either implicitly or 
explicitly, the reliability of the entire smooth potential approximation. Thus, 
e.g., it has been suggested [[l2| that "the approximation of a smooth potential 
is useful for studying orbits, but not for studying their divergence." This is of 
course a problematic statement in that the distinction between exponential and 
power law divergence, emblematic of the differences between regular and chaotic 
behaviour, lies at the heart of applications of nonlinear dynamics to galactic 
dynamics. If the Lyapunov exponents associated with the bulk potential have 
nothing to do with the TV-body problem, one must perforce reject completely all 
conventional applications of nonlinear dynamics to galactic astronomy. 

The crucial point, then, is that there does appear to be a well-defined contin- 
uum limit, even at the level of individual orbits |T^ , [r^ , |l5t . Suppose that a smooth 
density distribution, corresponding either to an integrable potential or to a po- 
tential admitting large measures of regular orbits, is sampled to generate a fixed, 
i.e. frozen in time, TV-body density distribution, and that the trajectories of test 
particles evolved in this frozen distribution are compared with smooth potential 
characteristics with the same initial conditions. In this case, there is a precise 
sense in which, as TV increases, the frozen-TV trajectories converge towards the 
smooth potential characterstics. Both visually and in terms of the complexity |l6| 
of their Fourier spectra, the frozen- TV trajectories come to more closely resemble 
the smooth potential characteristics; and, viewed mesoscopically, the frozen-TV 
and smooth potential orbits remain closer in phase space for progressively longer 
times. In particular, a frozen-TV orbit corresponding to an integrable character- 
istic will have a large Lyapunov exponent \n even if, visually, it is essentially 
indistinguishable from the regular characteristic! 

But how can this be? The key recognition here is that two 'types' of chaos 
can be present in the TV-body problem, characterised by two different sets of 
Lyapunov exponents associated with physics on different scales. Close encoun- 
ters between particles trigger microchaos, a generic feature of the TV-body prob- 
lem, which leads to large positive Lyapunov exponents xn- If, however, the 
bulk smooth potential is chaotic, one will also observe macrochaos, which is 
again characterised by positive, albeit typically much smaller, Lyapunov expo- 
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nents xs- Suppose, for example, that one compares the evolution of two nearby 
chaotic initial conditions in a single frozen-TV background or the same chaotic 
initial condition evolved in two different frozen- N realisations of the same bulk 
density. In this case, one typically observes a three-stage evolution, namely: (1) 
a rapid exponential divergence at a rate \n set by the true Lyapunov expo- 
nents associated with the iV-body problem, which persists until the separation 
becomes large compared with a typical interparticle spacing; followed by (2) 
a slower exponential divergence at a rate comparable to the (typically much 
smaller) smooth potential Lyapunov exponent Xs, which persists until the sep- 
aration becomes macroscopic; followed by (3) a power law divergence on a time 
scale ex (In N)to- For regular initial conditions, the second stage is absent and 
the time scale for the third stage scales instead as iV 1 / 2 ^. 

Microchaos becomes stronger as N increases in the sense that the value of 
Xn increases with increasing TV p7[ . Despite this, however, it becomes progres- 
sively less important macroscopically in that the range of the chaos, i.e., the 
scale on which the microchaos-driven exponential divergence of nearby orbits 
terminates, decreases with increasing N . In the limit N — > oo microchaos will 
become completely irrelevant but, for finite N, it does have an effect, at least on 
sufficiently short scales; and it is possible from an TV-body simulation to extract 
estimates of both xn and the typically much smaller xs p5| . 



3.2 Modeling discreteness effects as friction and noise 

It has been long recognised that, for sufficiently small N and/or over sufficiently 
long times, discreteness effects will not be completely negligible. Systems like 
galaxies are 'nearly collisionlcss' in the sense that the stars interact primarily 
via collective macroscopic forces associated with the bulk density distribution; 
but, at least in principle, if one waits long enough discreteness effects should 
have an appreciable effect. 

Astronomers are accustomed to modeling discreteness effects in the context of 
a Fokker-Planck description analogous to that formulated originally in the con- 
text of plasma physics |Tq] . However, it is not completely clear to what extent 
this is really justified. The conventional Fokker-Planck description was formu- 
lated originally to extract statistical properties of orbit ensembles and distribu- 
tion functions over long time scales, assuming implicitly that the bulk potential 
is regular. To what extent, then, can Langevin realisations of a Fokker-Planck 
equation yield reliable information about individual orbits over comparatively 
short time scales, particularly if the orbits are chaotic? 



Analyses of flows in frozen- N systems indicate 14 that, at the level of both 
orbit ensembles and individual orbits, discreteness effects can in fact be mod- 
eled extremely well by Gaussian white noise in the context of a Fokker-Planck 
description, allowing for a dimensionless diffusion constant D cx 1/N, consistent 
with the predicted scaling D cx InA/N, with A the so-called Coulomb loga- 
rithm For localised ensembles of initial conditions corresponding to both 
regular and chaotic orbits, phase mixing in frozen- N systems and phase mixing 
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in smooth potentials perturbed by Gaussian white noise yield virtually identi- 
cal behaviour, both in terms of the evolution of various phase space moments 
such as the emittance and the rate at which individual orbits in the ensemble ex- 
hibit nontrivial 'transitions', e.g., passing through some entropy barrier horn one 
phase space region to another. And similarly, a comparison of frozen- TV orbits 
and noisy smooth potential orbits with the same initial condition reveals that 
their Fourier spectra typically exhibit comparable complexities. Gaussian white 
noise is even successful in mimicking some of the effects of microchaos. If, e.g., 
one tracks the divergence of two noisy orbits with the same chaotic initial condi- 
tion evolved in a smooth potential, one observes the same three-stage evolution 
as for a pair of frozen- N orbits evolved in two different frozen- N potentials. 

An example of this agreement is illustrated in Fig. |2[ which exhibits data 
generated by averaging over 100 pairs of orbits evolved in frozen-A^ density 
distributions which correspond in the continuum limit to a triaxial homogeneous 
ellipsoid with axis ratios 1.95 : 1.50 : 1.05, perturbed by a spherically symmetric 
central mass spikes ('black hole'). The top two solid curves represent (from top to 
bottom) results for N = 10 4 5 and TV = 10 5 5 . The four dotted curves represent 
analogous results derived for pairs of noisy orbits evolved from the same initial 
conditions in the smooth potential with (from top to bottom) diffusion constant 
D = 10 -4 , 10~ 5 , 10 -6 , and 10 -7 . The near-coincidence of the top two solid 
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Fig. 2. The mean spatial separation between the same initial conditions evolved in 
two different frozen-iV backgrounds (solid curves) and different noisy orbits evolved 
in the smooth potential from the same initial condition (dots). The solid line has a 
slope 0.022, equal to the mean value of the smooth potential Lyapunov exponent xs- 
The dashed curve has a slope 0.75, equal to the mean value of the TV-body Lyapunov 
exponent XN- 
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and dotted curves indicates that discreteness effects for N = 10 p+1 / 2 are well- 
mimicked by Gaussian white noise with D = 10~ p . 

Such striking agreement suggests strongly that investigations of how orbits 
in smooth potentials are impacted by the introduction of friction and noise 
can provide important insights into the role of graininess in real galaxies. It is 
customary to assert that, in a system as large as a galaxy, discreteness effects 
reflecting close encounters between stars are unimportant because the relaxation 
time tn on which they can induce appreciable changes in quantities like the 
energy is orders of magnitude longer than the age of the Universe [[l9|. This is 
likely to be true if the galaxy is an exact equilibrium, especially an equilibrium 
characterised by an integrable potential. However, the assertion is suspect if (as 
must usually be the case) the system is only 'close to' an equilibrium or near- 
equilibrium, especially if the bulk potential is characterised by a phase space 
admitting a complex coexistence of regular and chaotic orbits. 

Over the past decade, analyses of flows in time-independent Hamiltonian sys- 
tems have revealed that even very weak perturbations, idealised as friction and 
white noise corresponding to tu ~ 10 6 — 10 9 tr> and, hence, D ~ 10~ 6 — 1CP 9 , 
can have significant effects within a time as short as lOOto or less by facilitating 
phase space diffusion through cantori or along the Arnold web [ ^o[|2l| , |22|| . The 
basic point is that the motions of chaotic orbits in a complex potential can be 
constrained significantly by topological obstructions like cantori or the Arnold 
web which, albeit not completely preventing motions from one phase space re- 
gion to another, serve as an entropy barrier to impede such motions. In many 
respects, the physical picture is similar to the elementary problem of effusion of 
gas through a tiny hole in a wall. There is nothing in principle to prevent a gas 
molecule from passing through the hole and, hence, escaping from the region 
to which it is originally confined; but, if the hole is very small, the time scale 
associated with this effusion can be extremely long. 

In the same sense, and for much the same reason, chaotic orbits trapped in 
one phase space region may, in the absence of perturbations, remain stuck in that 
region for a very long time. However, subjecting the orbits to noise will 'wiggle' 
them in such a fashion as to increase the rate at which they pass through the 
entropy barrier, thus accelerating phase space transport. Numerical simulations 
indicate that, in at least some cases, this escape process can be well approximated 
by a Poisson process, with the number of nonescapers decreasing exponentially 
at a rate A that is determined by the perturbation (2^,|2^]. This effect appears 
to result from a resonant coupling between the orbits and the noise. White noise 
is characterised by a flat power spectrum and, as such, will couple to more 
or less anything. If, however, the noise is made coloured, i.e., if instantaneous 
kicks are replaced by impulses of finite duration, the high frequency power is 
reduced; and, if the autocorrelation time becomes sufficiently long that there is 
little power at frequencies comparable to the orbital frequencies, the effect of the 
noise decreases significantly. Significantly, it appears that, overall, the details of 
the perturbation may be largely irrelevant: additive and multiplicative Gaussian 
noises tend to have comparable effects and the presence or absence of friction 
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does not seem to matter. All that appears to matter are the amplitude and the 
autocorrelation time upon which there is a relatively weak, roughly logarithmic, 
dependence. 

But what does all this imply for a real galaxy? Given that collisionless near- 
equilibria must be more common than true equilibria, it would seem quite possi- 
ble that, during the early stages of evolution, a galaxy might settle down towards 
a near-equilibrium, rather than a true equilibrium, e.g., involving what have been 
termed p5| 'partially mixed' building blocks. If discreteness effects and all other 
perturbative effects could be ignored, such a quasi-equilibrium might persist 
without exhibiting significant changes over the age of the Universe. If, however, 
one allows for discreteness effects or, alternatively, other perturbations reflect- 
ing, e.g., a high density cluster environment, the orbits could become shuffled in 
such a fashion as to trigger significant changes in the phase space density and, 
consequently, a systematic secular evolution p(J. 

Such a scenario could, for example, result in the destabilisation of a bar. Many 
models of bars {e.g. (^7)) incorporate 'sticky' |^8| chaotic orbits as part of the 
skeleton of structure, replacing crucial regular orbits which can be absent near 
corotation and other resonances. Making these 'sticky' orbits become unstuck 
could cause the bar to dissipate. Similar effects could also cause an originally 
nonaxisymmetric cusp to evolve towards a more nearly axisymmetric state. To 
the extent that the triaxial Dehnen potentials are representative, one can argue 
that chaotic orbits may be extremely common near the centers of early-type 
galaxies, but that many of these chaotic orbits are extremely sticky |^9| and, 
as such, could help support the nonaxisymmetric structure. Perturbations that 
make these sticky orbits wildly chaotic could de facto break the bones of the 
skeleton supporting the structure and trigger an evolution towards axisymmetry. 

4 Experimental tests of galactic dynamics 

4.1 Similarities between galaxies and nonneutral plasmas 

Even though electrostatics and Newtonian gravity both involve 1/r 2 forces, elec- 
tric neutrality implies that the physics of neutral plasmas is very different from 
the physics of self-gravitating systems. Viewed over time scales > t^, nonneutral 
plasmas and charged-particle beams are also very different from self-gravitating 
systems: the attractive character of gravity leads to phenomena like evaporation 
and core-collapse which cannot arise in a beam or a plasma. If, however, one re- 
stricts attention to comparatively short times ijj , much of the physics should 
be the same. Theoretical expectations, supported by numerical simulations, sug- 
gest that it is the existence of long range order, not the sign of the interaction, 
which is really important; but, to the extent that this be true, collisionless plas- 
mas and collisionless self-gravitating systems should be quite similar. 

Typical sources of charged-particle beams configure the beams in trains of 
'packets' or 'bunches', as they are termed by accelerator dynamicists. The ob- 
jective of a good high-intensity accelerator is to generate bunches comprised 
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of a large total number of charges confined to a small phase space volume and 
then accelerate those bunches to very high energies while minimising any growth 
in emittance. As one example, modern photocathode-based sources of electron 
beams routinely generate bunches comprised of some 10 10 — 10 11 electrons with 
transverse 'emittance' e of a few microns. (Here e = e/t>o, where vq is the mean 
axial velocity of the particle distribution.) The energy relaxation time i# asso- 
ciated with such bunches typically corresponds to the time required for a bunch 
to travel a distance ~ 1 km or so which, in many cases, is much longer than any 
distance of experimental interest, so that the bunches are 'nearly collisionless.' 

Models of equilibrium configurations of nonneutral plasmas and charged- 
particle beams confined by electromagnetic fields can be characterised by a 
complex phase space quite similar to that associated with models of ellipti- 
cal galaxies and, as such, have orbits with very similar properties. For exam- 
ple Q , the so-called 'thermal equilibrium model' Q of beam dynamics, which 
involves a self-interacting nonneutral plasma in thermal equilibrium confined by 
an anisotropic harmonic oscillator potential, is strikingly similar p9[| to the non- 
spherical generalisations of the Dehnen potential of galactic astronomy in terms 
of such properties as the degree of 'stickiness' manifested by chaotic orbits or 
how the relative measure of chaotic orbits and the size of the largest Lyapunov 
exponent vary with shape. 

As in galactic dynamics, questions have been raised regarding the validity of 
the continuum approximation for nearly collisionless charged particle beams ||3l| . 
However, comparatively short time integrations (t <C tu) involving discreteness 
effects and the nature of the continuum limit in nonneutral plasmas f^j yield re- 
sults essentially identical to what is observed for gravity - although the behaviour 
associated with neutral plasmas is very different. In particular, the macroscopic 
manifestations of phase mixing, for both regular and chaotic orbits, are indistin- 
guishable, and the coexistence of microchaos and macrochaos persists unabated. 

Nontrivial effects associated with a time-dependent potential have also been 
predicted for both nonneutral plasmas Q and charged particle beams g]. Al- 
though the time-dependence that is envisioned in a beam is typically less violent 
than that anticipated in violent relaxation within a galaxy, such is not always 
the case. Indeed, there is compelling experimental evidence that, in a beam, such 
a time-dependence can have the undesireable effect of ejecting particles from the 
core into an outerlying halo . 

Perhaps most interesting, however, is the fact that numerical simulations that 
reproduce successfully 'anomalous relaxation' observed in real laboratory exper- 
iments involving accelerator beams have shown compelling evidence of chaotic 
phase mixing. One classic example involves the propagation of five nonrelativis- 
tic high-intensity beamlets in a periodic solenoidal transport channel, where 
self-consistent space-charge forces are extremely important J35|. Ideally, these 
beamlets should exhibit coherent periodic oscillations (quite literally disappear- 
ing and reappearing) which might be expected to decay only on a relaxation time 
scale tji that corresponds to a propagation distance ~ 1 km. However, regard- 
less of how well the beam was matched to the transport channel, the beamlets 
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were seen to reappear only once, at a point ~lm from the source, disappearing 
completely within 2 m or so (c/. Fig. [6.10] in Their failure to reappear 

again would seem to reflect some collisionless process that, in effect, causes the 
particles to 'forget' their initial conditions. 

Detailed simulations using the particle- in-cell code WARP |3(|, which do an 
extremely good job of reproducing what is actually seen, demonstrate seemingly 
unambiguously that, because of the time-dependent space-charge potential, a 
large fraction of the particles in the beam experience the effects of strong, possi- 
ble transient, macrochaos |37| , ^8| . This is, e.g., evident from Fig. 3, which illus- 
trates the evolution of representative test particles which interact with the bulk 
potential but not with each other. Here the left hand panel shows snapshots of 
the beam after it has travelled distances m, 0.98 m, 2.88 m, 5.24 m, 11.52 m, 
and 31.68 m, with the representative ensembles superimposed. The right hand 
panel exhibits the evolution of the emittances e x and e y for these ensembles. It 
is evident that the initially localised ensembles are diverging exponentially so 
as to fill much of the accessible phase space, and that this exponential diver- 
gence coincides with the beamlets losing their individual identities. Also evident 
is the fact that the behaviour observed here is very similar to that exhibited in 
Fig. 1 which, recall, was generated for orbits in a perturbed Plummer potential 
exhibiting damped oscillations. 

4.2 Testing galaxy evolution with charged particle beams 

The aforementioned similarities between galactic astronomy and charged parti- 
cle beams suggest the possibility of using accelerators as a laboratory for astro- 
physics in which one can perform experimental tests of galactic dynamics, a pos- 
sibility currently being developed by a University of Florida - Fermilab/Northcrn 
Illinois University - University of Maryland collaboration. This collaboration, 
which has the dual aims of (1) obtaining an improved understanding of the ap- 
plicability of nonlinear dynamics to nearly collisionless systems interacting via 
long range forces and (2) using that understanding to generate more sharply 
focused bunches by minimising undesirable increases in emittance, is currently 
planning concrete experiments which can, and presumably will, be performed on 
the University of Maryland Electron Ring ( UMER) currently under construction. 
Here a number of obvious issues, all experimentally testable, come to mind: 

How ubiquitous is chaotic phase mixing as a source of anomalous relaxation? 
Older experiments with lower-intensity beams, where the space-charge forces 
were comparatively unimportant, tended not to manifest extreme examples of 
anomalous relaxation. Anomalous relaxation appears more common in high in- 
tensity beams, especially in settings where a time-dependent density distribution 
generates a strongly time-dependent potential; and it is obvious to ask whether 
chaos is the principal culprit. The idea here is to identify the types of sce- 
naria that tend generically to yield anomalous relaxation and to determine, e.g., 
whether such scenaria tend typically to be associated with a bulk potential that 
incorporates a strong, roughly oscillatory component. Do numerical simulations 
of orbits ensembles evolved in such systems exhibit evidence of chaotic phase 
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mixing? And do individual orbits in those ensembles exhibit strong exponential 
sensitivity, associated, e.g., with transient chaos? 

Do instabilities tend to trigger transient chaos? Instabilities in collisionless 
systems can exhibit behaviour qualitatively similar to that associated with tur- 
bulence in collision-dominated systems, but it is well known that turbulence is 
a strongly chaotic phenomenon. This possibility is especially interesting in that 
turbulence is another setting where different 'types' of chaos, characterised by 
wildly different time scales, can act on different length scales. 

What types of geometries, both strongly time-dependent and nearly time- 
independent, tend to yield the most efficient chaotic phase mixing and the largest 
measures of chaotic orbits? Do time-dependent evolutions involving strongly con- 
vulsive oscillations tend generically to exhibit especially fast relaxation? And do 
they tend to yield especially large amounts of chaos, as probed by the relative 
measure of chaotic orbits and/or the sizes of the largest (finite time) Lyapunov 
exponents? To the extent that bulk properties of such 'accelerator violent re- 
laxation' correlate with the degree of chaos exhibited in the evolving beam, and 
that the degree of chaos correlates with the form of the macroscopic time depen- 
dence, one will have a physically well-motivated explanation of which sorts of 
scenaria would be expected to exhibit complete and efficient violent relaxation 
and which would not! 

Is it, e.g., true that, for nearly axisymmetric configurations, prolate (or 
oblate) bunches tend generically to exhibit especially large amounts of chaos? 
And do any such trends that are observed coincide with trends observed in 
models of galactic equilibria [^9f ? Even if a beam bunch remains nearly axisym- 
metric during its evolution, the acceleration mechanism can - and in general will 
- change its shape as it passes down an accelerator in a fashion that depends 
on the accelerator design. The obvious question, then, is whether the oblate or 
prolate phase tends generically to be especially chaotic. 

Addressing this and related issues could provide important insights as to why 
galaxies have the detailed shapes that they do, a general question for which, at 
the present time, no compelling dynamical explanation exists. One knows, e.g., 
that elliptical galaxies tend to have isophotes that are slightly boxy or disky, 
and this boxiness or diskiness correlates with such properties as the rotation 
rate, the steepness of the central cusp, and the size of any deviations from ax- 
isymmetry [^] . Must all these effects be attributed to the detailed form of the 
formation scenario, or is there a clear dynamical explanation? Is it, e.g., true 
that the observed deviations from perfect ellipsoidal symmetry conspire to re- 
duce the relative number of chaotic orbits or to increase the numbers of certain 
regular orbit types required as a skeleton to support the observed structure? 

One might also use accelerator experiments to probe the role of discrete sub- 
structures and the extent to which they can be modeled as friction and noise in 
the context of a Fokker-Planck description j|0|. If the injection of a beam in- 
volves a large mismatch, a significant charge redistribution will occur, resulting 
in violent relaxation, 'turbulent' behaviour, and the formation of substructures 
('lumps') on a variety of scales. To the extent that such a time-dependent evo- 
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lution can be described in a continuum approximation, one might then expect 
that the bulk potential will correspond to a highly complex time-dependent 
phase space and that the substructures could act as a 'noisy' source of extrinsic 
diffusion, facilitating both transitions between 'sticky' and 'wildly chaotic' be- 
haviour and, in some cases, transitions between regularity and chaos. Given the 
evidence (c/. p3fl) that, at least over short times, discreteness effects act similarly 
for attractive and repulsive 1/r 2 forces, such insights could be directly related to 
such issues as the destabilisation of bars in spirals and/or the secular evolution 
of nonaxisymmetric ellipticals towards more nearly axisymmetric states. 

Do systems tend to evolve in such a fashion as to minimise the amount of 
chaos? There is an intuitive expectation amongst many galactic astronomers 
(c/. jll]]) that galaxies tend to evolve towards equilibria which incorporate few if 
any chaotic orbits, i.e., that nature somehow favors nearly-regular equilibria. It 
would certainly appear true that a model must incorporate significant numbers 
of regular orbits to support interesting structures like bars and/or triaxiality, 
but this does not a priori preclude the possibility of chaotic orbits also being 
present. Generic time-independent three-degree-of- freedom potentials are neither 
completely regular nor completely chaotic, admitting instead a complex coexis- 
tence of regular and chaotic phase space regions. The obvious question then is: 
are galactic equilibria or near-equilibria typically well-represented by potentials 
which are generic in this sense; or are they, for reasons unknown, special in that 
they tend to be rather nearly regular? 

5 Conclusions 

This paper has focused on several fundamental issues that arise in attempts to 
apply nonlinear dynamics to real galaxies, many-body systems characterised by 
a self-consistcntly determined bulk potential which, during their most interesting 
phases, can be strongly time-dependent. As recently as a decade ago these issues 
would have been considered of largely academic, rather than practical, interest. 
However, recent observational advances - which facilitate improved high resolu- 
tion photometry of individual objects as well as statistical analyses of large sam- 
ples with varying redshift - and improved computational resources - which allow 
unparalleled explorations of multi-scale structure -, together with the recogni- 
tion that the basic physics can be also probed in the context of charged particle 
beams, imply that theoretical predictions regarding these 'academic' issues can 
in fact be tested observationally, numerically, and experimentally. 

Thanks to Ioannis Sideris for providing Fig. 1 and to Rami Kishek for providing 
Fig. 3. Thanks also to Court Bohn for his comments on a preliminary draft of 
the manuscript. Limited financial support was provided by NSF AST-0070809. 
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Fig. 3. Evolution of five representative ensembles of test particles in the five-beamlet 
simulation. The left hand panel shows snapshots of the ensembles at (top-to-bottom 
left column) m, 0.98 m, 2.88 m, and (top-to-bottom right column) 5.24 m, 11.52 
m, and 31.68 m, with the x- and y-axes labeled in meters. The right panel shows the 
evolution of the logarithm of the emittances e x and e y as a function of distance S(z) 
along the accelerator. 



